The Hamiltonian boundary term and quasi-local energy flux 
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The Hamiltonian for a gravitating region includes a boundary term which determines not only the 
quasi-local values but also, via the boundary variation principle, the boundary conditions. Using our 
covariant Hamiltonian formalism, we found four particular quasi-local energy-momentum boundary 
term expressions; each corresponds to a physically distinct and geometrically clear boundary con- 
dition. Here, from a consideration of the asymptotics, we show how a fundamental Hamiltonian 
identity naturally leads to the associated quasi- local energy flux expressions. For electromagnetism 
one of the four is distinguished: the only one which is gauge invariant; it gives the familiar energy 
density and Poynting flux. For Einstein's general relativity two different boundary condition choices 
correspond to quasi-local expressions which asymptotically give the ADM energy, the Trautman- 
Bondi energy and, moreover, an associated energy flux (both outgoing and incoming). Again there 
is a distinguished expression: the one which is covariant. 
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PACS numbers: 04.20.Jb, 04.65. +e, 98.80.-k 

I. INTRODUCTION 

The Hamiltonian that generates the dynamical evolu- 
tion of any physical system within a (finite or infinite) 
spatial region necessarily includes, in addition to an inte- 
gral over the spatial volume, an integral over the bound- 
ary of the region. Our concern is with this Hamiltonian 
boundary term: we wish to understand both its role and 
proper form. 

The Hamiltonian is, of course, related to energy. Every 
physical system carries mass-energy and consequently in- 
evitably generates gravity, hence gravity plays the major 
role in our analysis. We consider only geometric grav- 
ity theories, and here more specifically only Einstein's 
general theory of relativity (GR). 

The essential necessity for, and the role of, the Hamil- 
tonian boundary term for asymptotically flat spaces in 
GR was first clearly discussed in a seminal work of 
Regge and Teitelboim They argued that in order 
for the Hamiltonian to be functionally differentiable on 
the phase space of asymptotically flat spatial metrics, 
certain boundary integrals over the 2-sphere at spatial 
infinity must be included. These turned out to be in fact 
just the ADM asymptotic expressions for the con- 
served total quantities: energy, momentum and angular 
momentum (plus an additional expression for the center- 
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of-mass). Subsequently certain improvements were made 
in the analysis by Beig and 6 Murchadha Q; more re- 
cently Szabados has made some further refinements 

Although the total conserved quantities are well de- 
fined for asymptotically flat spaces, as is well known, 
there are no well defined local densities for these quanti- 
ties. This can be understood in terms of the equivalence 
principle, which precludes the detection of gravity at a 
point (see, e.g. p §20.4). The modern idea is that one 
should have quasi- local quantities: i.e., quantities asso- 
ciated with a closed 2-surface (for a nice review of this 
topic see 0). 

For finite regions, we found that the value of the 
Hamiltonian boundary term can determine the quasi- 
local quantities. Using essentially the same principle 
as used by Regge- Teitelboim (this principle will play 
the major role in our subsequent discussion) we found 
that the Hamiltonian boundary term also determines the 
boundary conditions. Using this idea along with our co- 
variant Hamiltonian formalism, for each dynamic field 
we found four types of quasi-local energy-momentum 
boundary term expressions; each corresponds to a physi- 
cally distinct and geometrically clear boundary condition 
Z 8, 9]. 

However a consideration of the phase space asymp- 
totics necessary for a well defined Hamiltonian naturally 
leads to the recognition that the conditions are not met 
in the radiating regime. This apparent difficulty provides 
the opportunity for our main new result. We show here 
how a fundamental Hamiltonian identity naturally leads, 
for each Hamiltonian boundary expression, to an associ- 
ated quasi-local energy flux expression. 
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We note that gravitational energy flux (aka balance) 
expressions have a long history, going back to Trautman 
[Tol |. Bondi [13 and Sachs Nevertheless there still 
is considerable interest in this topic and some significant 
progress has been made recently; see in particular ,,13i . il4 
[15L \Wi UtI IT^ . We include here some applications of 
our new natural energy flux expressions. 

An application of our formalism with its quasi-local 
energy and energy flux expressions to the electromagnetic 
field, where we can readily interpret the expressions and 
boundary conditions, is included to illustrate the ideas 
in a familiar setting. For the electromagnetic case one of 
our four expressions is distinguished: the only one which 
is gauge invariant; it gives the usual energy density and 
Poynting flux. 

When applied to Einstein's general relativity two dif- 
ferent boundary condition choices correspond to quasi- 
local expressions which asymptotically give the ADM en- 
ergy, the Trautman-Bondi energy and, moreover, an asso- 
ciated energy flux (both outgoing and incoming). How- 
ever once again there is a distinguished expression: in 
this case it is the one which is covariant. 

The plan of this work is as follows: In section II we 
discuss the first order Lagrangian formalism. In section 
III we consider local translational invariance. The basic 
Hamiltonian formulation is considered in section IV. In 
section V we consider refined boundary terms. There fol- 
lows a discussion of the phase space and the asymptotics 
in Section VI. Our new flux expressions are presented 
in section VII. In section VIII we consider the applica- 
tion of these ideas to electromagnetism. The application 
to Einstein's gravity theory is covered in section IX. A 
discussion forms the concluding section. 



Its variation has the form 

SC = d(6ip Ap)+6ipA — + — ASp, (2) 
dip op 

implicitly defining a pair of first order Euler-Lagrange 
expressions, which are explicitly given by 
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Sp 



-i,dp - d^A, 



(3) 



where <r := (—I)'''. The integral of C associates an action 
with any spacetime region. The variation of this action 
is given by the integral of SC. The total differential term 
in ^ then leads to an integral over the boundary of 
the region. Hamilton's principle — that the action should 
be extreme with Sip vanishing on the boundary — yields 
the field equations: the vanishing of the Euler-Lagrange 
expressions Q. 



III. LOCAL TRANSLATION INVARIANCE 

The action should not depend on the particular way 
points are labeled. Thus it should be invariant under dif- 
feomorphisms, including infinitesimal diffeomorphisms — 
which correspond to a displacement along some vector 
field N. From a gauge theory perspective such displace- 
ments are regarded as a "local translation". Under a 
local translation quantities change according to the Lie 
derivative. Hence, for a diffeomorphism invariant action 
the relation ((SJ should be identically satisfied when the 
variation operator S is replaced by the Lie derivative £n 
(= diN + i^d on the components of form fields): 



II. THE FIRST ORDER LAGRANGIAN 
FORMALISM 

Our formalism is intended to be applicable to general 
types of fields. Technically we find it convenient to repre- 
sent dynamic fields in terms of differential forms. (More 
precisely to accommodate all the fields found in nature we 
would use a collection of tensor and spinor valued forms 
including Dirac spinors and gauge potential one-forms 
along with the spacetime orthonormal coframe and con- 
nection one-forms). Here we develop the representative 
case of an k-form field ip (the field may take its values in 
some vector space and may thus carry some indices which 
are here suppressed; the generalization to include several 
fields, possibly of different grades, is straightforward). 

We proceed from the action principle. It can be shown 
that any action principle can be rewritten in an equiva- 
lent form, which (following e.g. ^^J) we refer to as first 
order, this is the most convenient form for our purposes. 
A first order Lagrangian 4-form for a k-form field ip has 
the form 



SC SC 

diiqC = £nC = d{£NpAp) + £Nf^^ — \'^^£nP- (4) 

Sp Sp 

This simply means that C is a 4-form which depends on 
position only through the fields </?, p. (For this to be the 
case the set of fields in C necessarily includes dynamic 
geometric variables, which means gravity.) 
From it directly follows that the 3-form 



niN) £Nip/\p-iNC 



satisfies the identity 



S C SC 

dUiN) EE £np a — — a £nP, 
Sp Sp 



(5) 



(6) 



showing that it is a conserved "current" on shell (mean- 
ing: when the field equations are satisfied). Substitut- 
ing into lO leads to the explicit expression H{N) = 
dliNP A p) -I- ^iNf A dp + ^dp} A iNP + in A, from which 
one can see that this conserved Noether translation cur- 
rent can be written as a 3-form linear in the displacement 
vector plus a total differential: 



C = dp Ap — A{p,p). 



(1) 



n{N) =: N^'H^,+dB{N). 



(7) 



3 



Compare the difFerential of this expression, dH{N) = 
dNf' AHf, + Nf^dHf,, with ® ; equating the dA^^ coeffi- 
cient on both sides reveals that 



S£ 5C 
rtf, = -VV' ^ + ''"^ ^ VP- 



IS) 



This identity is a necessary consequence of local diffeo- 
morphism invariance (i.e., symmetry for non-constant 
N'-'-). From this relation one can see that vanishes 
on shell; hence the value of the conserved quantity asso- 
ciated with a local displacement N and a spatial region S 
is determined by a 2-surface integral over the boundary 
of the region: 



E{N,^):= [ n{N)= (f B{N). 

JT. JOT, 



(9) 



The value is quasi-local. For any choice of N this expres- 
sion defines a conserved quasi-local quantity. 

What do these values mean? In general we do 
not have a clear physical interpretation. However, at 
least if the geometry on the boundary is not so far 
from flat space, for a suitable timelike(spacelike) quasi- 
translation displacement on the boundary the expression 
defines a quasi-local energy (momentum), and for a suit- 
able quasi-rotation(boost) it defines a quasi-local angu- 
lar momentum(centcr-of-mass). (Here we do not explore 
the important question of how to specifically make these 
quasi-displacement choices for a general region in order 
to obtain good physical quasi-local values.) 



IV. THE HAMILTONIAN FORMULATION 

From the first order field equations Q, by contrac- 
tion with a "time evolution vector field" N and using 
ZArrf</5 = £n^ — diN'P, iNdp = £nP — diNP, we get a 
pair of Hamiltonian-like evolution equations for the "time 
derivatives": £-M^p, £nP- A key identity involving these 
time derivatives is revealed by comparing two relations; 
on the one hand take the projection of 5L Q along N: 



ij^SC = iNd{5if /\p) + In 



5ip A 



S(p Sp 



A Sp 



£n{S(p Ap) ~ diiq{5ip Ap) 



-In 



Sip A 



5C_ 5C 

6(p Sp 



A Sp 



(10) 



on the other hand take the projection of the Lagrangian 
4-form iNjC, which from ^ is just £N'f Ap — H{N), and 
vary: 

SiNC EE S{£nV^p) - Sn{N) (11) 
EE S{£n<p) Ap+ £Nf ASp- 5n{N) 
= £nS(p Ap+ £n'P ASp - SHiN) 
= £N{Sip Ap) - Sip A £nP+ £nV ASp- Sn{N). 



Since is not varied the two relations are identical: 
SinC = inSC; consequently. 



Sn{N) = -Sip A £nP+ £nP ASp 
+diN{Sf Ap) — In 



(12) 



sc sc , 

Sip A- — h A dp 
Sip Sp 



The last term vanishes "on shell" . This relation identifies 
the Noether translational current Ti{N) as the Hamilto- 
nian 3-form (i.e., density), as the following considerations 
show. The integral of T-C{N) over a 3-dimensional region. 



ff(7V,S) := J^nN), 



(13) 



is the Hamiltonian which displaces this region along N, 
since the integral of its variation: 



sh{n,j:) = [ sniN), 



(14) 



yields, from (|13|) . (on shell) the Hamilton equations: 



£n'p 



sn{N) 

Sp 



£np ■ 



5U{N) 

Sip 



(15) 



if the boundary term in the variation of the Hamiltonian 
vanishes. In this case that means when Sip vanishes on 
dYi. Technically the variational derivatives of the Hamil- 
tonian H{N, E) are only defined for variations satisfying 
this boundary condition. In other words this Hamilto- 
nian is functionally differentiable on the phase space of 
fields satisfying this particular boundary condition. 



V. REFINED BOUNDARY TERMS 

In some important cases the fields of physical inter- 
est do not satisfy the aforementioned boundary condi- 
tion naturally inherited from the Lagrangian (the main 
example is the spacetime metric for an asymptotically 
flat region). A suitably modified formulation is needed 
to deal with this situation. One alternative is to mod- 
ify the Lagrangian 4-form itself by a total differential. 
As discussed in some detail in , this would modify the 
boundary condition on the whole 3-dimensional bound- 
ary of the spacetime region, thus inducing the same type 
of modification on the spatial boundary at large spatial 
distances as on the initial time hypersurface. However 
we actually want the freedom to adjust the boundary 
condition on the 2-dimensional boundary of the space- 
like region 9S] independently of the type of initial con- 
ditions imposed within S. For this purpose we focus on 
the Hamiltonian. We note that the Hamiltonian ||7J) has 
two distinct parts; each plays a distinct role. The proper 
density iVTi^, although it has vanishing value on shell, 
generates the equations of motion, whereas the boundary 
term B{N) determines both the quasi-local value (O and 
the boundary condition. Now the boundary term can be 
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adjusted — without changing the Hamilton equations or 
the conservation property © — indeed we can replace the 
2-form B{N) — i^ip Ap inherited from the Lagrangian 
by any other. Such an adjustment is in one respect just 
a special case of the conserved Noether current ambi- 
guity (i.e. for any 2-form Xi J and J' := J + dx are 
both conserved currents (3-forms) if dJ — 0, even though 
they define different conserved values). However here, in 
this Hamiltonian case, any such adjustment modifies — 
in parallel — not only the value of the quasi-local quanti- 
ties but also the spatial boundary conditions. Thus the 
boundary term ambiguity is under physical control: each 
distinct choice of Hamiltonian boundary quasi-local ex- 
pression is associated with a physically distinct boundary 
condition We thus find that the Hamilto- 

nian density always takes the general form 

Hb.cXN) ^ Nf^n^ + dBB.cXN) . (16) 

The subscript "B.C." here refers to the particular choice 
of built in boundary condition. 

In order to accommodate suitable boundary conditions 
we found that, in general, one must introduce (at the 
minimum actually only on the boundary, but to simplify 
the discussion here we presume it to be on any desired 
neighborhood of the boundary) certain reference values 
p, <f, which represent the ground state of the field — the 
"vacuum" (or background field) values (this is neces- 
sary in particular for fields whose natural ground state 
is non- vanishing, e.g. the spacetime metric). We take 
our boundary terms to be linear in Aip :^ ip — (p and 
Ap := p — p, so that they (and thus all the quasi- local 
quantities) vanish if the fields take on the ground state 
(i.e. reference) values. We presume that the reference 
values (like A'') are not varied: S(p = and Sp = 0, conse- 
quently 5A(p = Slp, 5Ap = Sp. Note that we presume the 
reference values (as well as N) to be defined on the dy- 
namic spacetime in the region of interest, independently 
of the dynamic fields; they are regarded as being fixed 
prior to (and thus independently of) the choice of E and 

5 = as. 

From our investigations 0,13 '^^ found two "covariant" 
boundary term alternates to B{N) := i^ip Ap, namely 

'BDirichict(A^) := ijv<P A Ap - ^A<y9 A ZArp, (17) 

eumaiin 

(N) :— ijsup A Ap — (,Atf AinP- (18) 

A short calculation (of the form STLb.c. — 5{T-L — 
dB)+SdBB.c.) shows that the variation of the associated 
Hamiltonian 3-forms (|16|1 have (on shell) the respective 
forms 

^WDirichlet(A^) = -S(p A £nP + £nV 

+diNiSip AAp), (19) 

SH-NcumanniN) = Sip A £ nP + £ N V ^ 

-diNiAip ASp), (20) 

corresponding to holding fixed certain (after integration 
certain projected) covariant sets of components, respec- 



tively the value of the field and the value of its canonically 
conjugate momentum; whence our labels. 

Moreover, we found two more physical interesting 
choices Q 

Sdynamic(^) := iN<f /\ Ap - qAlfi A tNp, (21) 
Sconstraincd(A^) := tNV /\ Ap ~ qAip A iNP- (22) 

The variation of the associated Hamiltonian 3-forms have 
(on shell) the indicated forms 

^Wdynamic(^) = Sip A £ nP + £n^ A Sp (23) 

+d{':,Sip A inAp — ipfAip A Sp), 

^'Hconstraincd(^) = Sip A £ NP + £n^ A Sp (24) 

+d{SiNip A Ap — qAp A Si^p), 

corresponding to holding fixed certain (after integration 
certain projected) components of respectively the "dy- 
namic parts" (the spatial puUback) of ip, p and the "con- 
strained parts" (the "time" projections) ijv<P, iNP- 

In each of these cases the boundary term in the Hamil- 
tonian variation has a certain symplectic structure which 
pairs certain control and response quantities . Within 
each of our two sets of expressions the pairs are simply 
related by an interchange of "control" and "response", 
formally 5 ^ A, A — > —S. 

Note that one expression stands out: for fields which 
allow trivial reference values, (p = = p, the bound- 
ary term i3dynamic(-^) vauishcs. These fields, with this 
choice of boundary condition, make no explicit contribu- 
tion to the quasi-local boundary term. Thus there is a 
certain preferred boundary expression — and thus a pre- 
ferred boundary condition — for this large class of fields, 
a class which includes all the necessary physical fields — 
aside from dynamic geometry gravity. 

An instructive discussion of boundary conditions as- 
sociated with variational principles appears in Lanczos 
j23 |. section 11.15. The variational principle always gives 
us the correct number of boundary conditions. However 
it should be remarked that it cannot guarantee that the 
boundary conditions are the proper ones for the exis- 
tence and uniqueness of solutions; that would depend on 
the particular quality and type of the equations, which 
is influenced especially by the metric signature (our for- 
malism "formally" does not care) and the details of the 
Lagrangian "potential" K{ip,p). Our general formalism 
cannot take such particulars into account. 

Our philosophy is that normally one should "control" 
on the boundary the indicated variations. Lanczos also 
discusses an interesting alternative: one could instead 
take "free" variations on the boundary. To have the 
boundary term in the variation vanish with free varia- 
tions of the "control" variables one must then require 
that the associated "response" vanishes; this yields what 
is referred to as the natural boundary condition. For our 
expressions it would simply mean that the response fields 
would take, on the boundary, the pre-specified "refer- 
ence" values. Hence, as far as imposing boundary con- 
ditions is concerned, one can achieve the same result by 
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these two different approaches: (i) control the variable 
to the desired boundary value, or (ii) take free variations 
with the appropriate reference fields chosen to have the 
desired boundary values. Note, however, that there are 
differences in the resultant quasi-local values. In partic- 
ular a free variation for TiDiiichict requires from H19|l that 
p = p. Consequently the value of the Hamiltonian is de- 
termined by a reduced form of the boundary expression 
(|17|l : —'^AipAiNP- Similarly, free variation for TiNcumann 
leads to (f = if and the reduced boundary term i^^tpAAp. 
On the other hand, for T^dynamicj free variations lead to 
the vanishing of i^Atp and ii^Ap, which yields no formal 
reduction of /Bdynamic, while free variation for Ticonstrainod 
leads to Ap = = A</j, consequently ^constrained van- 
ishes! Hence we find a curious fact: for free variations 
there is boundary expression which is distinguished by 
having a vanishing value, it is the expression which is 
complimentary to the one which vanishes for trivial ref- 
erence values. 

Lanczos has discussed, via a detailed example, how one 
could do more complicated things, where certain vari- 
ables are "controlled" and others are varied "freely" . In 
all cases the boundary variation principle yields the cor- 
rect number of boundary conditions. For our expressions 
there are many possible combinations of "free" and "con- 
trolled" variations. Using the "natural boundary condi- 
tions" associated with "free variation" are an interesting 
option which merit further exploration. Based on our 
present understanding, however, for most applications we 
favor purely "controlled" variations. 



VI. THE PHASE SPACE AND ASYMPTOTICS 

For finite regions (which is our primary interest) these 
boundary terms in the variation of the Hamiltonian tell 
us exactly what needs to be held fixed (not constant but 
rather "controlled" , i.e. the function on the boundary is 
predetermined by some outside agent). For asymptoti- 
cally flat regions, however, it is not sufficient to just say 
we want the field to vanish at infinity. Rather one should 
take into account the asymptotic fall off rates. Asymp- 
totically we want to allow the variations and responses 
to be like the differences between generic solutions. The 
various boundary terms we have constructed will all en- 
able the Hamiltonian to be well defined on the phase 
space of fields with asymptotic behavior for all typical 
physical fields. Detailed investigations have been done in 
particular for Einstein's gravity theory, general relativity 
(GR) beginning with the pioneering work of Regge and 
Teitelboim P|. This work was later improved by Beig 
and 6 Murchadha Q and more recently further refined 
by Szabados j^. Here we give a simplified summary of 
certain relevant conclusions of these works. 

It turns out that one should impose different fall off 
rates for the terms with different parity. For the fields 
it is sufficient to take the respective asymptotic fall offs 



and parities to be 

Aip^O++ , Ap « + 0+ , (25) 

where r'^^O^j can be finite as r — > oo and ± indicates the 
parity. Here, for all types of tensors and forms, we define 
the parity to mean the parity of the components in an 
asymptotically Cartesian reference frame dx"^ . Now the 
Cartesian components of the 2-surface area element have 
odd parity, consequently even parity 2-forms automati- 
cally have vanishing 2-surface integral. 

For asymptotically flat spaces the displacement should 
asymptotically be a Killing vector, i.e. an infinitesimal 
Poincare displacement. More precisely one can allow 

iV^ « (ao + a^)^ + (Ao + A+)^x^ (26) 

where ctg is a constant translation parameter and Ag'' — 
a||"^' is a constant asymptotic infinitesimal Lorentz 
boost/rotation parameter; the even parity part of the 
perturbations in these parameters can have the indicated 
l/r fall off, but any odd parity perturbation to A should 
fall off faster. Note that the A^ correction, because of its 
coordinate coefficient, means that it can be regarded as 
an odd parity perturbation of a. Thus we can, without 
any change in the conserved total values or vanishing of 
the boundary term in the Hamiltonian variation, allow 
super-translations — i.e., asymptotically non-constant (= 
angular dependent) terms of finite magnitude — as long 
as they are odd parity. Even parity supertranslations, 
on the other hand, would change the value of the quasi- 
local quantities and would also yield in general a non- 
vanishing boundary term in the Hamiltonian variation. 
If one really wanted them one could allow even parity 
supertranslations, but only at the expense of requiring 
boundary conditions on the fields more strict than H25|) . 

With the asymptotics H25I26|I it is straightforward to 
verify, as we just stated, that all four of our Hamiltonians 
are differentiable on the specified phase space — since all 
the boundary terms in the variations of the Hamiltonians 
vanish asymptotically; moreover the quasi-local expres- 
sions all give the same (since they differ by terms of the 
form NAipAp which vanish asymptotically) finite con- 
stant values (independent of the perturbation af and 
super-translation Af x) for the energy-momentum and 
angular momcntum/center-of-mass — even for fields, like 
the metric or frame, with « 1. 

Note that the indicated asymptotic behavior is suffi- 
cient but hardly necessary. In practice only the specific 
projected component combinations that actually show up 
in the boundary integrand need be so restricted (more- 
over one need only require this behavior up to a closed 
form). However it is not possible to give a general co- 
variant formula for such details; for such refinements one 
must examine how each component for a particular field 
in a specific theory actually occurs in the expression. 

We also note that one could even turn things around 
and take the finiteness of the quasi-local expressions and 
the vanishing of the boundary term in the variation of the 
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Hamiltonian as defining a norm that would determine the 
largest possible phase space. We do not pursue this idea 
here. 

At spatial infinity the aforementioned asymptotics are 
physically reasonable. Let us now consider what can 
be expected if the boundary of our 2-surface d'E ap- 
proaches null infinity. One can imagine it following the 
characteristic propagation surfaces of outgoing radiation. 
Long range radiation fields (e.g. electromagnetism) have 
slower fall offs, like Ap k, df ~ Oi. Then it may seem 
that we will have a serious problem, since the bound- 
ary terms in the variation of our various Hamiltonians 
will not vanish, so the Hamiltonian is no longer func- 
tionally differentiable. This seeming calamity is actually 
providential — as was recognized long ago (see e.g. the 
remark on page 160 in 125(|) but not investigated until 
more recently — it is directing us to additional physics 
contained within the formalism, namely energy flux ex- 
pressions. 



VII. FLUX EXPRESSIONS 

Given any vector field M = d/dX one could calculate 
(on shell) the associated change in any quasi-local quan- 
tity directly from the boundary expression: 

■^Eb.cXN,^) = I £mHb.c.{N) = i £mBb.c.{N). 
dX Jy. Jos 

(27) 

In this fashion one could calculate from such a "fiux" ex- 
pression, for example, the change in the quasi-local linear 
momentum under a rotation or the time rate of change 
of the angular momentum. In particular this approach 
can be specialized to calculate the "time rate of change" 
of the "energy" associated with the "time displacement" 
N = d/dt itself 



£NnB.C.{N) 



as 



£nBb.c.{N). 



(28) 

The 2-surface integral then defines an "energy flux" . 
(Note: we use the labels "time" and "energy" here for 
convenience in our description, however the actual mean- 
ing is the change along N in the conserved quantity as- 
sociated with N] e.g. we could consider a rotation: N = 
d/d(j)). Of course there is another natural way to calcu- 
late the "time rate of change" of any conserved quantity, 
namely one simply rearranges the conserved current for- 
mula, such as d^j'^ = to Stj" = —dkj^, and then in- 
tegrates over a spatial region; the spatial divergence via 
the divergence theorem yields a boundary 2-surface flux 
integral. Such expressions are often referred to as "bal- 
ance equations" . Expressed in terms of a closed 3-form 
they follow just from integrating £Nj = diMj- In par- 
ticular for energy this leads (on shell) to £n'Hb.c.{N) = 
diNnB.cXN) = diNdBB.cXN) = d£NBB.cXN)- Thus 
this alternative approach also leads to H28|) . 



For the fiux of "energy" (and only for the flux of energy, 
i.e. the change of E{N) along N , not for the "time rate of 
change" of any other quantity like linear or angular mo- 
mentum, etc.) there is another way of calculating — the 
analog of the classical mechanics calculation (for conser- 
vative Hamiltonian systems) of 



5H = (fSpk-pkSq^ 



E:= H = (29) 



under the replacement S d/dt, where the remarkable 
cancellation is a consequence of the particular (symplec- 
tic) form of the Hamilton equations. Specializing the 
relations H19I2UI24I25|) with S ^ £n (we are presuming 
that £N<f, £nP vanish, that could be a strong restric- 
tion on the reference or on N), we note that the 3-form 
parts vanish identically (again this is associated with the 
symplectic form of the Hamilton equations), hence the 
respective boundary fiux expressions are the integrals of 

£N'Hl)irichlct{N)=diNi£NV A (30) 
£N'HNcumann{N)=diN{-Aip A £nP), (31) 
£N'HdymimiciN)^d{q£Nip/\iNA.p-iNA(pA£Np),i32) 
£AfHconstraint(A^)=d(i7v£Af<;5AAp-^A<y9AZAri;7vp).(33) 

These expressions hold even if does not vanish (i.e, 
when we have global but not local translation symmetry). 
Note that they are significantly different in appearance 
from those obtained using if^ along with H17I18I21I22|I : 

i^ArHDirichlct(A^) d£N{iN^ A Ap - <:,Aip A ZArp),(34) 
£ NT-iyienmnnniN) = d£N{iN'fi A Ap - i;A(fi A iArp),(35) 
"£AfHdynamic(A^) = d£N{iN'f A Ap - ^Aifi A iArp),(36) 
£ NTi-constrmntiN) = d£N{iNV ^ Ap - qAlf A ZArp),(37) 

which hold only as long as Ti^ vanishes. Via a straight- 
forward but non-trivial calculation using the field equa- 
tions H15() it can, of course, be explicitly verified that the 
two respective forms are actually equivalent — when 
vanishes. 



VIII. APPLICATION TO 
ELECTROMAGNETISM 

To illustrate these ideas first in a familiar setting let us 
examine (vacuum) electromagnetism. One could consider 
this as a source field for gravity, but it is more instructive, 
and still sufficient for our needs here, to more simply just 
consider vacuum electromagnetism in Minkowski space. 
In that case the formalism developed above, with the 
important exception of the "on shell" vanishing of 7i^, 
is still applicable. The first order Lagrangian 4-form for 
the (source free) U(l) gauge field one-form A is 



C^M^dAhH- A H, 

ZAq 

yielding the pair of first order equations 

1 



dH = 0, 



dA 



A, 



-*H = 0. 



(38) 



(39) 
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These are just the vacuum Maxwell equations with *H = 
XqF XodA; hence H = -Xo*F, and d*F = 0. (Here 
Ag"^ is the vacuum impedance.) With N = dt and the de- 
composition A = (—(f),Ak) we find that inF = incLA — 
£ nA—cHnA corresponds to Fok = Ak+dk<p = —Ek- The 
magnetic field strength is Fij :— diAj — djAi —: EijkB'^. 
Hence Hoi — —Xa^Foi = —XoBi, Hij — _~'Xa*Fij _ = 
—XoeijkE^ . The natural reference choice is A = = H. 
The Hamiltonian 3-form is 

^Bx;.(^) = -iNAdH - cLAMnH 
1 



+iN[^''HhHj+dBB.c.- (40) 

In conventional tensor index notation the volume density 
part is 



n 



EM 



<j)duE^ + ]^{d,A,~djAi)F'= 



--E'^Ek - -F'^F,,, 



1 

4' 



(41) 



After eliminating the 2nd class constraint, Fij = diAj — 
djAi, it corresponds to the familiar Ao[^(i?^ + B^) — 
(f)dkE^]] the scalar potential acts as a Lagrange multiplier 
to enforce the Gauss constraint dkE^ = 0. 

For the four considered boundary choices we have 



jEM 
-"Dir 
rEM 
"Ncu 

EM 



InAH = \a<t)E''dSk, (42) 
A A iNH = -XoAiBje'^'^dSk, (43) 



B. 



0, 



dyn 

Bf^ = iNAH + AAinH 

= Xo{cf>E'' ~ A,Bje'^'')dSk. 



(44) 



(45) 



The Hamiltonian variations have the respective forms 

Sn%^{N) = field equation terms 

+d{SiNAH -SAMnH), (46) 
'^'^Neu(^) = Reld equation terms 

+d{-iNA5H + AASiNH), (47) 
STi^^iN) = field equation terms 

+d{-iNASH -SAAinH), (48) 
6T-lf^{N) — field equation terms 

+d{diNAH + AASiNH). (49) 

Here our interest is not in the field equations but in 
the total differential term which, upon integration, be- 
comes a boundary term indicating the boundary condi- 
tion. Briefly: the choice corresponds to fixing the 
scalar potential and the components of the vector poten- 
tial parallel to the 2-surface (the gauge independent part 
of the latter fixes the normal component B± of the mag- 
netic field). H™, enables the fixing of certain projected 
components of F, namely E± and a part of _B||. The 



choice Hfj^ is used for fixing E± and B±, while is 
used for fixing the scalar potential and a part of B^. 

The physical meaning of such boundary conditions are 
well known especially in the electrostatics case. Fixing 
E± on the boundary corresponds to fixing the surface 
charge density. Here is an instructive physical applica- 
tion: use a battery to first charge up a capacitor which 
contains a dielectric which can be inserted/removed, then 
disconnect the battery and measure the work needed to 
remove/insert the dielectric (in the process the potential 
will vary but the charge is fixed, no current or power will 
flow into/out of the system, the system is adiabatically 
insulated). Alternately leave the battery connected and 
measure the work needed to displace the dielectric — then 
the potential is fixed although the charge will now vary, 
so in this latter case current and hence power flows into 
or out of the system. The respective boundary terms 
in the variation of the Hamiltonian are —(f)5E^dSk and 
5(f)E^dSk- The point we wish to emphasize is that both 
boundary condition choices are physically meaningful; 
they correspond to real situations actually encountered 
in practice. Nevertheless there is a preferred choice. 

One expression stands out: the choice is the only 
one in which the value of the Hamiltonian is gauge in- 
variant. Moreover, in addition to the neat property of en- 
joying a vanishing boundary term, it has an extra virtue 
of considerable physical importance: namely the Hamil- 
tonian density is non-negative. Consequently the associ- 
ated energy has a lower bound and the system has a nat- 
ural vacuum or ground state: zero energy for vanishing 
fields. The value of the Hamiltonian can be inter- 

preted as the internal energy, whereas the other expres- 
sions include some additional energy on the boundary 
of the system associated with maintaining the boundary 
condition. 

The respective energy flux expressions from (|30H33f) 
are 

£N'H^ = d{iNdiNAAH-diNAAiNH~iNFAiNH) 
= Xod [{d^AoBj-^dtAoe^jkE'' -E,Bj)dx' Adx'l{50) 

£N'H§fu = d{AAiNdiNH -iNAAdiNH) 
^ Xod [iAQ^^BJ - A.dtBj) dx' A dx^] , (51) 

£nH^^, = -d{iNF AiNH) 



Xod {-E,Bj dx' A dx') , (52) 

^ d{iNdiNAAH + AAiNdiNH) 
Xod [- {^dtAoe.^kE'' -f A,dtB,) dx' A dx^ . (53) 



(Note we cannot calculate the correct energy flux in this 
case by expressions like H28l34llH7|l unless we include an 
additional non- vanishing contribution from Ti.^.) All of 
these Hamiltonians and their associated energies are re- 
ally describing the same physical laws (the differences 
in the right hand sides just correspond to differences in 
the left hand sides), however all except are bound- 

ary condition choices which are gauge dependent descrip- 
tions. Clearly the choice, associated with fixing 
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the normal components E± and B± on the boundary, is 
preferred; it is the one suitable for most physical appli- 
cations. It gives the usual energy density and Poynting 
flux. 



IX. APPLICATION TO EINSTEIN'S GRAVITY 
THEORY 

Einstein's gravity theory, general relativity (GR), can 
be formulated in several ways. For our purposes the 
most convenient is to take the orthonormal coframe 
= i}>^i^dx^ and the connection one-form coefficients 
uj^p — V"' pi^dx^ as our geometric potentials. Moreover 
we take the connection to be a priori metric compati- 
ble: Dgap := dgap —uj'^a9-yp —^'^pga'j = 0. Restricted to 
orthonormal frames where the metric coefficients are con- 
stant, this condition reduces to the algebraic constraint 

We consider the vacuum (source free) case for simplic- 
ity (the inclusion of sources is straightforward). GR can 
be obtained from the first order Lagrangian 4-form 

Cgr = ^''''APc.p+D^>'AT^-V''''A{po^l3-^r,^l3), (54) 

where f2"^ := dw"^ -t- w"^ Ato^ p is the curvature 2-form, 
D'd'^ := dd^ + to'^^, A i?" is the torsion 2-form, and we 
have made use of the convenient dual form basis rj"'!^--- ■= 
*(z9" Ai}l^ A ■■■). The 2-forms fl"''^, V"''^ and p^o are 
antisymmetric. We take k := SttG/c"^ = Stt. 

In this Lagrangian the frame and connection conjugate 
momenta t^, pafj and the auxiliary field V"^ all function 
like Lagrange multiplier fields; their variation imposes, 
respectively, the torsion vanishing condition Di}^ ~ 0, 
the multiplier value V"^ = fi"^, and the conjugate mo- 
mentum value Pa/3 = (2K)~^r]af3. The connection one- 
form variation gives (in vacuum) 

Dpap + i5[/3 A T„] = 0. (55) 

Since Dpap oc Drjais = D'd'' AT]ai3^ = we get i^ipATa] = 
0, from which it follows that = (in vacuum). The 
frame variation gives 

Dr^ + ^y"^ ^Vaf^^ = Q- (56) 

Substituting the already determined values for and 
yafS yields fl"^ A rj^p^ = 0, the vanishing of the Einstein 
tensor 3-form, i.e. the vacuum Einstein equation. 

For most physical fields we can get by with a trivial ref- 
erence. However in the case of GR we certainly need the 
refinements of both adjusting the Hamiltonian boundary 
term "by hand" — as was first clearly argued by Regge and 
Teitelboim (RT) [l| — and introducing a reference. With 
just the boundary term in the Hamiltonian naturally 
inherited from the Lagrangian, B{N) = ijsfU)"^ A p^fB, 
the boundary term in the variation of the Hamiltonian 
has the form i^iSuj"^ A rjap), which does not vanish for 



asymptotically flat fall offs. A simple reason for introduc- 
ing the reference values, at least for the connection, is to 
render covariant the manifestly non-covariant "natural" 
Hamiltonian boundary term. (By the way the need for 
reference values for GR was not apparent in the ADM 
and RT formulations; there it was hidden in the choice 
of asymptotically Cartesian coordinates. The explicit 
need for a reference metric in the asymptotic Hamilto- 
nian boundary term was, to our knowledge, first clearly 
apparent in the work of Beig and 6 Murchadha Isj). 

In the vacuum case (or more generally as long as our 
boundary is in the vacuum region) the frame conjugate 
momentum field = 0, so only the connection and 
its conjugate momentum make explicit contributions to 
our gravitational quasi- local expressions. (When sources 
are included we can always choose the boundary condi- 
tions so that the source fields have vanishing quasi-local 
boundary term, however via the gravitational field equa- 
tions the sources, of course, indirectly influence the values 
of the gravitational field variables on the boundary and 
thereby do contribute to the quasi-local values.) 

Now one of our dynamic gravitational variables is not 
a tensor field, so there are certain terms in our quasi- 
local expressions which include the non-covariant factors 
i^uj"'!^ or i^uj"''^. As discussed in more detail in 0,0], 
the physical contribution due to these terms is obtained 
by replacing them by dI^^A^"! or l)['^A^"l. Taking these 
considerations into account, along with the identification 
Pap = (2K)~^7]a/3, we obtain our quasi-local GR Hamilto- 
nian boundary term expressions for gravitating systems 
(in vacuum regions): 

^™(^) := 7^ [^u^'^^AiNriap + D^^N^^^^Tlap] ,(57) 

BT{N) ■■= ^ i^uj'^^AiNfiap+D^PN^^^^riap] ,(58) 

S^y^(iV) := ^ i^^uj'^f'AiNfiap + D^^N^^^^riap) ,(59) 

Cn(^) := 7^ (^uj^^AiNria0 + D^I'N°^^/\ijap] .(60) 

As already mentioned these quasi-local expressions will 
certainly yield finite values when integrated over an 
asymptotic 2-sphere for asymptotic Killing displacements 
of the form (|25|l . at least when the variables satisfy the 
asymptotic parity and fall off conditions 

(5t?«0++02", 5lj f=i O:^ + Ot ■ (61) 

(At null infinity some parts of the connection actually 
have slower fall off yet, as we shall see, the quasi-local 
integrals are all still finite.) 

The variations of the associated Hamiltonians have the 
respective forms 

Sn'i^iN) = field equation terms (62) 

+ ^d{-iNl^uj'''^ 5r)ap + Aw"^ A 5iNr]ap): 
Zk 

6n™iN) = field equation terms (63) 
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Ik, 



field equation terms 



(64) 



6T-lf^{N) = field equation terms 



(65) 



Here our interest is not in the field equations but in 
the total differential term which, upon integration, be- 
comes a boundary term indicating the respective bound- 
ary conditions. Briefly: Ti.^^ requires fixing (after inte- 
gration: certain projected components of) the orthonor- 
mal coframe, while requires fixing (certain projected 
components of) the connection. Whereas Hdyn associ- 
ated with fixing the spatial projections of the frame and 
connection, and 7i^^ is associated with fixing the time 
components of the frame and connection. The bound- 
ary terms in these Hamiltonian variations vanish asymp- 
totically (spatially) with the aforementioned fall off and 
parity conditions. Asymptotically at spatial infinity the 
quasi-local quantities obtained from all four of our ex- 
pressions are compatible with the analysis of Beig and 
O Murchadha or more precisely with the refinement 
thereof of Szabados . 

The associated energy flux expressions, calculated ac- 
cording to the respective prescriptions H3UH33|I . (presum- 
ing that £ffd = = £n'^, i-S- -/V is a Killing field of the 
reference geometry) take the form 

(66) 

£Nn2'^ iN) = ^d(£Arijvw"'^ Ar?a/3 - £ nco'^'^ AinAtJo^p), 

(67) 

£Nn^^^iN) = ^di^iNAio''^ £Nria.p-£Nio'"^ AlNAj^cti), 



(68) 



£Nn'^^^{N)^—d{£NiNLo^''Ar^^0 + ALo''PA£NiNV<.p). 
Zk, 

(69) 

Only one, £]\[H'^^ (|66|l . is free from non-covariant fac- 
tors. 

Of course one should check the actual values given by 
these expressions. We have reexpressed (|57I66|I in the 
asymptotically null regime using the NP spin coefficient 
formalism and obtained good results ,26.] . Moreover we 
are presently working on adapting our expressions to con- 
formal infinity; the results of that investigation will be re- 
ported in due course. Meanwhile, similar to the calcula- 
tions in ^ , we have tested all of these energy and energy 
flux expressions at null infinity on the full Bondi-Sachs 
[l^ form of the metric using Reduce and EXCALC. In 
particular, for the value of the expression (which had 
been calculated earlier |23) we obtain results essentially 



identical to those reported in Jl6| obtained using SHEEP 
to evaluate the Freud holonomic expression (which is the 
expression Trautman used in his original work 0|). This 
is not at all surprising, since this particular orthonormal 
frame expression, for any frame satisfying the asymptotic 
"no rotation" gauge condition "diak] ~ O2, is asymptot- 
ically equivalent to both the Freud superpotential 0, 
and the expression of Katz and coworkers 14]; the latter 
also gives good results at null infinity. 

Our three other boundary expressions give similar but 
not identical results. Both for brevity and to more clearly 
show the main similarities and differences we here present 
only the essential details for the simpler, axi-symmetric 
Bondi metric ■ We take the coframe to be of the form 



^* = e'^+'^'du + e^-'I'dr, 



r = e'^-'I'dr, 

,9" =re^{de -r-^Udu), d"^ = re"'^ sin6ld^, (70) 



where e^"^ = 1 - 2m{u,9)/r + O2, 7 = c{u,9)/r + O2, 
j3 = -c^/Ar"^ -h O4, and [/ = -dec - 2c cot 9 + d. For 
the reference we take (j) = ^ = (3 = U = Q. For N = du 
we find for the asymptotically contributing part of the 
energy expressions 
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= B^^^ = — [4m sin 9 + 2de{U sin 9)] d9 A dip, (71) 
2k 



B2^ = B2^ = —[4:{m + ccu) sm9+2d0{Usm9)]d9 Adip. 
" 2k 

(72) 

Upon integration over the sphere the terms with 
96i(?7sin0) make a vanishing contribution due to the reg- 
ularity conditions on c, d^c at the poles (note: similar 
terms appear in the 2-surface integrands in [T^ 12^'). 
Then the two expressions in H71|l give the standard Bondi 
mass aspect, while the other two do not. 

The asymptotically contributing part of the associated 
quasi-local flux values were found to be 

£Qjif^ = £a„H°^ = ^d \-^cl^M^9d9 A dp\ ,(73) 

^djij" - ^dJilfr. = [4cc„„ sinfl d0 A d^\ .(74) 



Thus two expressions (|75|) directly give the standard 
Bondi flux loss. The other two expressions are actually 
describing exactly the same physics albeit for a different 
energy expression, as can be easily seen from comparing 
the time derivative of the expression H71|) equated to 
with the time derivative of the expression H72|l equated to 
((Tljl . Note that one can equally well compute an incom- 
ing flux from past null infinity by assuming a dependence 
on the advanced time w ~ i + r in place of the retarded 
time u ~ t — r. 

From these calculations we can see that the functions 
/3 and IJ do not play a major role in this limit. This ob- 
servation justifies the following simpler calculation, which 
still includes enough of the essential qualities while show- 
ing that our expressions capture both the incoming and 
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outgoing quasi-local flux. Let us evaluate the energy 
and energy flux expressions for the following simplified 
asymptotic form of the orthonormal coframe: 

§v ^ e~'<rsm9dip, (75) 

where = (j){t,r) — 0{l/r), 7 — 'y{t,r) — 0{l/r). For 
the reference we take = = 7. The connection one- 
form components are 

^tr ^ (jj'e^'l'dt- ^e-^'I'dr, 

= rje~'l'+'fd9, 

uj^'^ = -r7e~"*~'^sin6'dv3, 

cj^^ = -{l + r-i')e-'l'+^de, 

uj^f = co&edip. (76) 

For our expressions with N — dt the DN terms make no 
asymptotic contribution. The key factors are 

Aw'''^ = [l-(l + r7')e'^+^] dO, 

t^uf'P = \V-(\-ri)e'^^'^\widd^. {11) 

From these we find the asymptotically contributing part 
of the quasi-local energy boundary expressions: 

Sf^ = = -!-4to sin 6'd6' A d^, (78) 

= -!-4[to - r27(7„ - 7„)] sin dO A dip, (79) 
2k 

where, as usual, m := (r/2)(l — e^"^) and u := t — r, 
V := t + r. Upon integration over the 2-sphere r = const 
we find that two boundary conditions, corresponding to 
the boundary terms B'^^ and B^^, give the Bondi mass. 
The associated quasi-local fiux values are found to be 

£nH^^ = £Nnf^ = [Ar^jj'smOde A dip] 

= [4r2 (72 - 72) sin 9 dO A d^] , (80) 

i^nf" = £Nn^^, = ^d [-Ar^jj' sinede a d^] 

= i-d [-4^2^(92^ _ dl-f) sine d9 A dip] .(81) 

Integrating (|78I - I81|I over a large 2-sphere at constant r, 
we find that all control modes are consistent with a Bondi 
news type energy flux loss and gain from both outgoing 
and incoming radiation: 

m = r^J' = r^{-fl-jl). (82) 

However only the Hamiltonians and give this 
relation directly. For most applications "H^^ would be 
the preferred choice, because the associated flux relation 
(|66|l is free of non-covariant factors. 



X. DISCUSSION 

The Hamiltonian for a gravitating system in a flnite or 
infinite region necessarily includes a boundary term. We 
have inquired into the significance and best form of this 
boundary term. We found that it determines not only 
the quasi-local values but also, via the boundary varia- 
tion principle, the boundary conditions necessary for a 
well-defined Hamiltonian. Wc noted that it is always 
possible (and necessary, at least for gravity) to include 
in it non-dynamic reference values for the dynamic vari- 
ables. Using our covariant Hamiltonian formalism and 
an identity associated with the variation of the Hamil- 
tonian, we have identified, for each dynamic field, four 
special quasi-local energy-momentum boundary term ex- 
pressions; each corresponds to a physically distinct and 
geometrically clear boundary condition. We showed how 
a fundamental Hamiltonian variation identity naturally 
forces us, for radiating asymptotics, to relax the well- 
defined Hamiltonian requirement and thereby obtained 
the associated quasi-local energy flux expressions. When 
the formalism is applied to electromagnetism one of the 
four is distinguished by gauge invariance (it gives the 
familiar energy density and Poynting flux). When the 
formalism is applied to Einstein's general relativity, two 
different boundary condition choices correspond to quasi- 
local expressions which asymptotically give the ADM en- 
ergy, the Trautman-Bondi energy and, moreover, an as- 
sociated energy flux (both outgoing and incoming), but 
once again there is a distinguished expression: in this 
case it is the one which is covariant. 

Here we make a few further remarks. First we men- 
tion that our ideas regarding variational principles, sym- 
plectic structure and the role of boundary terms have 
been much influenced by several sources, especially Tul- 
czyjew and Kijowski [2^|2|| and his coworkers. Although 
we considered radiation using the Hamiltonian, there are 
many aspects associated with the Hamiltonian in the ra- 
diating regime that we did not discuss; these issues are 
nicely treated in The uniqueness of the Trautman- 
Bondi mass is established in [ig . 

We, like many others, took a Hamiltonian approach. 
While some of our expressions are similar to those that 
appear in other approaches, e.g. the Noether charge ap- 
proach, our formalism includes certain features that are 
unusual. Features which distinguish our formalism in- 
clude: We have made extensive use of differential forms — 
because of their convenience for integration over domains 
with boundaries, the space-time decomposition, as well 
as the representation of geometric and gauge fields. We 
started from a first order Lagrangian — mainly because it 
facilitates the passage from the Lagrangian to the Hamil- 
tonian but also because of its linearity re coupling via 
connections. We used the (co)frame and connection as 
independent variables. The Hamiltonian formalism we 
have developed is 4-covariant. 

Our aim has been to find 4-D covariant 2-form and 3- 
form expressions which can be integrated over any suit- 
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able desired regions; essentially we want to find a good 
Hamiltonian density for the region of interest. In this 
approach the particular 3-surface E and its boundary 
d'E are quite incidental. Accordingly, our reference val- 
ues and displacement are independent of these surfaces. 
Note also that we do not decompose our expressions into 
intrinsic and extrinsic parts with respect to such sur- 
faces, as was done both traditionally and in many modern 
works, e.g. |3^ wherein some of the themes consid- 
ered here have been treated with the aid of such decom- 
positions. 

Aside from some important clarifications of the for- 
malism, the main new ingredients here involve our vari- 
ational identity, the relaxation of the boundary varia- 
tional principle, and the flux expressions associated with 
our four quasi-local expressions, along with the fact that 
in each application one expression is distinguished (by 
being gauge invariant for the relevant type of gauge 
covariance — U{1) or Lorentz in our electromagnetic and 
gravitational examples, respectively). Although in our 
discussion we referred to energy and energy-flux, what 
we really considered is the value of the Hamiltonian for 
a prescribed displacement N; the associated flux expres- 
sion that gives the rate of change of this value along N. 
Our formal results may have other useful applications 
besides energy. 

Localization of energy-momentum has been an out- 
standing problem from the very beginning of GR. Tradi- 
tional methods, both Noether spacetime translation sym- 
metry and decomposition of the field equations, have only 
led to a variety of reference frame dependent expressions 
(generally referred to as pseudotensors) for the energy- 
momentum density. Thus in addition to the ambiguity of 
which expression should be preferred, there was also the 
ambiguity of the choice of reference frame. The bound- 
ary variation principle along with the Hamiltonian for- 
malism tamed these inherent ambiguities, giving them a 
clear physical and geometric interpretation: the choice of 
expression is related to the choice of boundary conditions 
and the choice of reference frame is associated with the 
choice of ground state ,9, 2L, 2^ . Within the covariant 
Hamiltonian formalism using this principle we had iden- 
tified certain special "covariant symplectic" Hamiltonian 



boundary term quasi-local expressions 0i IM 0- Here, 
based on a fundamental Hamiltonian identity, we have 
identified the associated energy flux expressions. More- 
over we found that, for both source fields and gravity, 
among the four one particular quasi-local Hamiltonian 
boundary term and its associated flux expression was dis- 
tinguished. 

An additional virtue enjoyed by the respective distin- 
guished boundary terms H44I57() is that the Hamiltonian 
for this quasi-local choice has positive energy, not only 
for elcctromagnetism (as we discussed above), but also 
for gravity. More precisely (for acceptable choices of 
the variables) the gravitational Hamiltonian with the dis- 
tinguished boundary term (|57|l is non-negative, at least 
when S is maximal or nearly maximal. Here we just 
make mention of two lines of argument that lead to this 
conclusion. One can obtain this result by adapting to fi- 
nite regions the global positivity proof ^31] using the SOF 
gauge conditions ,32] ; one can instead use an adaptation 
of the Shi and Tam proof , which guarantees the pos- 
itivity (for mean convex 2-surfaces) of the Brown- York 
[s^ l quasi- local expression (the latter, as shown in j^, 
agrees with for certain choices). We are presently 
working on a detailed account of these two arguments. 
Of course it would be nice to have a stronger and more 
general energy-momentum positivity proof for our dis- 
tinguished covariant-symplectic quasi-local GR Hamilto- 
nian boundary term. 
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